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SECTION  1 


INTRODUCTION  AND  SUMMARY 


1.1  PURPOSE 

The  purpose  of  this  study  is  to  develop  a  multi-sensor,  multi-target  data  association  technique 
applicable  to  ASW  platforms  utilizing  long  range  passive /active  sonars  in  a  convergence  zone  ocean 
environment. 

1.2  GENERAL  PROBLEM  DESCRIPTION 

Advanced  passive  sonar  systems  can  detect  and  track  targets  at  long  range  while  operating  in  a 
convergence  zone  (CZ)  multi-target  environment.  These  tracks  are  not  continuous,  hence  inter /intra 
CZ  fading  gives  rise  to  the  need  for  Temporal  Association.  Since  these  tracks  may  be  viewed 
simultaneously  by  different  sensors,  Mutual  Association  is  also  required. 

For  many  operations,  only  passive  sonars  are  used;  however,  provision  must  be  made  for  active 
sonar,  surface  radar  and  sensors  on  other  /supporting  platforms.  With  passive  operation,  bearing  and 
frequency  tracks  are  assumed  to  be  smoothed  to  obtain  at  least  the  mean  bearing/frequency  and 
corresponding  rates. 

Two  target  types  are  common;  submarines  and  surface  ships.  Typical,  North  Atlantic,  surface 
shipping  tracks  of  12-hour  duration  are  shown  in  Figure  1-1.  A  typical  bearing-time  history  for  a  long 
range  passive  sonar  in  the  circled  area  is  shown  in  Figure  1-2.  The  gaps/fades  in  the  tracks  are  due  to 
the  high  propagation  loss  that  exists  between  convergence  zones.  The  long-term  data  association 
problem  is  to  bridge  these  gaps  using  the  bearing-time  data  shown  as  well  as  other  available  measurables 
such  as  bearing  rate,  frequency  rate,  active  sonar  range/bearing  and  surface  (radar)  data. 

1.3  BACKGROUND 

In  1977,  a  GE  study  sponsored  by  NOSC  t  *  1  addressed  the  data  association  problem  for  a  surface 

ship  employing  hull  and  towed  arrays.  The  development  of  the  Expected  Likelihood  approach  started 

at  that  time.  A  two-stage  look  back  approach  was  simulated  using  the  data  of  Figure  1-2.  Although 

many  assumptions/simplifications  were  made,  the  approach  appeared  to  be  better  matched  to  the 

problem  than  other  known/existing  techniques.  Essentially  all  of  the  existing  techniques  assume  good 

observability  and  Maximum  Likelihood  related  measures.  SI  Chou  discusses  the  Expected  Likelihood 

(31 

approach  in  a  recent  NOSC  report 
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Concurrently  with  the  above  study  ^  l  a  data  association  process  was  developed  in  the  GE 

simulation  facility  for  surface  ship  ASW.  It  was  based  on  a  weighted  sum  of  geometric,  frequency  and 

classification  correlation  indices  (not  based  on  likelihood  techniques).  Recently  some  of  the  bearing, 

frequency  techniques  described  in  this  report  were  compared  with  this  process  and  found  to  be  an 

[2] 

improvement 

1.4  SUMMARY 

Starting  with  probability  of  association  as  the  fundamental  association  measure  (utilizing  prior 
target  and  environment  models),  the  Expected  Likelihood  Ratio  emerges  as  the  key  measure.  This 
measure  is  the  weighted  average  of  the  likelihood  function  with  averaging  taking  place  over  elements 
of  the  target  state  vector;  the  weights  are  the  prior  probability  density  function  of  these  elements. 

Complexity  was  reduced  by  careful  choice  of  the  state  variables  and  appropriate  simplifications. 
For  long  term  (CZ  to  CZ)  association,  the  best  variable  choice  was  considered  to  be  range/bearing  at 
two  time  points  (reference  time  (tQ)  and  candidate  time  (tj)).  The  averaging  process  then  reduced  to  an 
integration  over  only  two  variables;  the  range  pair  (R0Rj).  For  short  term  or  simultaneous/mutual 
association  the  preferred  state  variable  set  was  range,  bearing,  target  speed  and  heading.  The  reduction 
again  resulted  in  a  two  variable  integration  process;  range  and  heading. 

Thus,  two  Expected  likelihood  Ratio  algorithms  evolved:  (1)  long  term,  applicable  when  At  is 
long  (e.g.,  greater  than  30  minutes)  and  (2)  short  term  (e.g..  At  less  than  30  minutes). 

The  long  term  (CZ  -  CZ)  algorithm  is  given  in  Section  4.  In  the  numerator  note  that  array 
gain/propagation  loss  is  factored  into  the  detection  term,  target  speed/heading  priors  in  the  target  term, 
geometric  measurements  in  the  Fy  term  and  frequency  spectra  in  the  last  term.  Integration  is  over  the 
range  pair  (RQRj)  where  RQ  is  range  at  the  time  of  the  reference  segment  (o)  and  Rj  is  range  at  the  time 
of  the  jth  candidate.  The  denominator  is  a  relatively  simple  normalizer  with  one  integral  for  each 
segment.  The  short  term  algorithm  is  similar  in  form,  but  differs  in  detail  as  described  in  Appendix  A6. 

1.5  ASSUMPTIONS 

a)  All  targets  are  of  the  same  general  class  and  hence  no  distinction  is  made  between  surface 
ship  and  submarine  targets. 

b)  Targets  are  non-maneuvering. 
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c)  Measurement  errors  are  gaussian  and  independent. 

d)  A  priori  target/environment  statistics  are  known. 

e)  Geometric  and  frequency  spectrum  data  are  the  only  measurables/attributes;  classification, 
demon  data,  etc.,  are  not  included. 

f)  Segment  estimates  of  mean  bearing,  bearing  rate,  etc.,  used  as  inputs  to  the  association  process 
have  previously  been  correctly  correlated /associated  and  bad  data  edited  out. 

1.6  CONCLUSIONS  AND  RECOMMENDATIONS 

The  Expected  Likelihood  technique  developed  in  this  study  provides  a  practical  approach  to  the 
problem  of  associating  data  segments  with  diverse  attributes;  geometric  and  frequency  spectrum  were 
included  but  others  can  be  added.  Prior  target  statistics  (e.g.,  speed/heading)  and  detection/environmental 
prediction  (e.g.,  probability  of  detection  versus  range)  are  incorporated  with  few  restrictions  on  the  form 
of  the  statistical  models.  The  formulation  does  not  depend  on  having  a  range  estimate  and  other  initial 
conditions  and  hence  will  work  long  before  localization  solutions  become  available.  The  formulation  is 
general  enough  to  include  the  mutual  and  temporal  association  of  data  segments  from  different  platforms 
and  hence  applicable  to  a  wide  variety  of  sonar /radar  sensors  including  single/multiple  ship  hull/towed 
arrays,  offboard  sonobouy /arrays  as  well  as  radar  and  other  surface  sensors. 

Recommendations  for  further  study  are  as  follows: 

a)  For  a  given  current  segment  a  likelihood  ratio  is  computed  for  each  past  segment,  one  at  a 
time;  a  so  called  single  stage  process.  A  significant  improvement  is  expected  when  the  two 
stage  (two  past  segments)  process  is  implemented. 

b)  Additional  attributes  should  be  included  in  the  likelihood  ratio  formulation  including  at 
least:  Demon  signatures  and  classification  decisions.  The  current  study  assumed  that  all 
targets  are  of  a  given  class,  namely  surface  ships  at  constant  course  and  speed;  maneuvering 
targets/subs  should  also  be  included  in  the  process. 

c)  Another  important  improvement  is  to  determine  how  to  make  better  association  decisions 
given  the  likelihood  ratio  measures.  Some  work  was  done  in  the  current  study  but  more  is 
requ.  d. 
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d)  Trade  studies  should  be  conducted  so  that  design  decisions  can  be  made.  A  good  example  of 
this  is:  should  the  amplitude  or  clipped  frequency  spectrum  be  implemented?  This  is  a 
performance  vs  complexity  trade. 

e)  The  ASW  simulator  should  be  further  exercised  to  provide  evaluations  of  the  association 
techniques  for  trade  studies,  technique  tuning  and  assessing  the  payoffs  for  improvements. 

f)  A  study  should  be  performed  to  determine  the  information/data  rate  currently  available 
from  existing  sensors/platforms/data  links.  The  payoff  for  more  data  and/or  additional 
sensors  should  be  determined  to  provide  planning/requirements  inputs. 
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SECTION  2 


FORMULATION  OF  ASSOCIATION  MEASURES 

As  a  target  is  tracked  through  a  convergence  zone,  the  resulting  smoothed  data  set  is  defined  as  a 
segment.  For  a  typical  ith  segment  at  time  tj,  the  data  is  defined  by  Z-  or  Zj(t-).  Three  types  of  data 
are  included:  Detection  (Dj),  Geometric  { Yj)  and  Frequency  Spectrum  (Sj)  so  that 

VID^]  (1.1.) 

A  current  or  reference  segment  is  identified  by  subscript  o  (i.e.,  ZQ),  a  candidate  segment  set  by 
subscript  j  (i.e.,  Zj)  and  the  remaining  segments  by  subscript  k  (i.e.,  Zk).  Thus,  the  total  data  set  is 
defined  by  the  alternative  forms. 

Z=  Zojk=  fZojZkl  *  lZo  Zj  Zkl 

In  general,  the  candidate  segment  set  can  include  more  than  one  segment.  These  definitions  are 
illustrated  in  Figure  2-1. 

The  problem  is  to  find  the  best  match  between  a  current /reference  segment  and  various  candidate 
segment  sets  (j  =  1,2,...). 


•EARING 
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Consider  the  general  association  problem  of  determining  “how  well”  a  current  data  segment  (ZQ) 
matches  with  the  jth  candidate  (Zj).  The  measure  of  association  is  given  by 


Foj  ■  P<VZ> 


(2.1) 


This  is  the  probability  of  association  of  the  o  and  jth  segment  set  conditioned  on  observing  the 
complete  segment  set  (Z). 


As  shown  in  Al^),  application  of  Bayes’  rule  and  manipulation  yields 


Apj  P  (AQj) 

Foj  =  p(Z)  /  P(Z  |  A) 


(2.2) 


The  key  element  is  the  Expected  Likelihood  Ratio  (A0j)  and  hence  the  focus  of  much  of  the 
subsequent  development.  The  other  numerator  term  P(AQj)  is  the  a  priori  probability  of  association 
of  the  current/candidate  set  (ZQj).  The  denominator  is  a  function  of  the  total  data  set  (Z)  and  hence 
is  a  normalizing  constant  during  subsequent  comparison/thresholding  of  the  candidate  sets.  P(Z)  is 
the  unconditional  density  function  of  Z.  A  is  the  condition  that  none  of  the  segments  are  associated. 


If  the  formulation  allows  for  only  a  single  past  segment  then,  as  developed  in  A1 


(2.3) 


where  Pg  is  the  probability  that  a  past  detection  occurs  and  m  is  the  number  of  past  segments  observed. 


(1)  Throughout  this  report  Al,  A2,  Bl,  etc.  are  appendices.  For  instance,  A2  is  Appendix  A, 
paragraph  2.  When  an  appendix  equation  is  referenced,  it  is  designated  by  (  ),  i.e.,  A2  (4). 
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LIKELIHOOD  RATIO 


The  expected  likelihood  ratio  for  the  current  o  and  jth  segment  set  is  given  by  (refer  to  Al) 

p(z0ji  M 

A°j  =  P(ZojIAoj) 

To  simplify  notation  the  oj  subscripts  will  be  omitted  and  n,  d  used  to  designate  the  numerator 
denominator  so  that 


(2.4) 


A  =  An/Ad  (2.4a) 

This  expresses  the  a  posteriori  probability  of  receiving  ZQ-  with  the  condition  that  the  segments  of 
ZQj  have  a  common  origin  and  hence  are  associated  (A0j).  The  non  association  condition  ( Aqj )  applies 
to  the  denominator. 


A  state  vector  (X)  is  used  to  predict  the  measurables  (Z). 


The  approach  used  to  develop  the  numerator  and  denominator  is  to  write  the  joint  measurement/ 
state  probability  density  function  (pdf)  and  then  integrate  out  the  state  variables,  thus  as  shown  in  A2, 


A 


V(2°‘ 


I  X)  P(X)  dX 


[fP(Zo  I  XQ)  P(XQ)  dXQ]  •  lfP(Zi  I  Xj)  P(Xj)  dXj] 
XQ  Xj 


(2.5) 


where  the  association  condition  is  implied  for  all  P(  ). 


The  conditional  density  functions  (i.e.,  P(ZQ  j  XQ))  are  likelihood  functions  and  the  others 
(i.e.,  P(XQ))  are  the  a  priori  (prior)  probability  density  functions  of  the  state  vectors  (i.e.,  XQ). 

P(ZQj  I  X)  is  the  product  of  the  likelihood  functions  of  each  measurable  in  the  oj  segments.  The 
denominator  integrands  apply  to  each  segment.  If  the  jth  candidate  consists  of  multiple  segments  then 
the  Xj  integral  is  replaced  by  the  product  of  multiple  integrals;  one  for  each  segment  in  the  candidate 
set. 
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NUMERATOR  -  An 

Consider  first  the  numerator  of  (2.5).  The  state  vector  (X)  consists  of  Geometric  and  Spectral 
(frequency)  elements  and  hence  will  be  expressed  as  X  =  [Xq,  Xg  J .  The  geometric  state  variables 
selected  for  long  term  (CZ-CZ)  association  are  given  by 

XG  =  IVoRj*jl  (2.6) 

where, 

RO0O  =  Range  and  bearing  at  reference/current  time  (tQ) 

RjO j  =  Range  and  bearing  that  applies  at  the  time  of  the  most  recent  segment  of 

the  jth  candidate  set  (tj). 

Using  the  form  of  Z  given  by  (1.1)  the  numerator  of  the  likelihood  ratio  as  developed  in  A3  is 
given  by 


(2.7) 


where  Roj0oj  are  shorthand  notation  for  RqRj  and  respectively,  primed  state  vector  notation 
(X'G)  means  that  the  state  variables  0Q  and  have  been  replaced  by  the  measured  values  (<T0<Tj). 


The  terms  of  (2.7)  are: 

•  p(XG)  —  A  priori  pdf  of  target  state  (Xq). 

•  P^Doj*XG^  —  Joint  Probability  of  detection  of  each  segment  conditioned  on  Xq.  This 

factors  in  the  environment  (propagation  loss/convergence  zones). 

•  Fy  —  The  geometric  measurement  likelihood  function  given  RQRj. 

•  pf®oj^oj  XG)  —  The  frequency  spectrum  likelihood  function  given  Xq. 


Note  that  all  of  the  non  gaussian  terms  are  outside  the  0Q0j  integral  thereby  allowing  an  analytic 
evaluation  of  Fy.  This  is  valid  for  the  long  term  (CZ-CZ)  association  but  not  for  short  term  association 
(e.g.,  within  a  single  CZ). 
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For  short  term  including  mutual  association,  refer  to  A6. 


DENOMINATOR  -  Ad 

Now,  consider  the  denominator  of  (2.5)  designated  by  A  d .  It  is  the  product  of  terms  A  .  where 

ai 

the  form  of  each  term  is  identical  except  for  the  subscript  identifying  the  segment  (i=o  j). 


Consider  a  typical  oth  segment  having  measurables  06  f.  As  developed  in  A4, 


Ad  =  P(0o)  J P(R0)  P(D0IR0)  P(0ofoIDoRo)  dR0]  P(S0ID0)  (2.8) 

°  Ro 

where  p (0Q)~  A  priori  bearing  state  pdf  evaluated  at  the  measured  bearing. 

P(0ofoIDoRo>  =  pdf  of  bearing  rate  and  frequency  rate  conditional  on  detection  and  range. 

Since  Of  (predicted)  are  a  function  of  cross  range  rate  (X),  in  addition  to  R, 
it  is  computed  by  (o  subscript  is  implied) 


P(0f]DR)  =  jf P(0f|DRX)  P(X)  dX 
X 


(2.9) 


COMPLETE  A 


Thus,  An  is  given  by  (2.7)  and  Ad  obtained  as  the  product  of  terms  like  (2.8).  Substitution  into 
(2.5)  yields  the  complete  expected  likelihood  ratio  given  by  (IT  is  product  notation) 


J  /  P(X’q)  P(D0Dj|Xq)  Fy  Aq(Xq)  dR0  dRj 
Ro  Rj _ 

n  /  PfRjflj)  P(Djl R|)  Pt^jfjlDjRj)  dRj 
i  Rj 


where  the  expected  likelihood  ratio  of  the  frequency  spectrum  data  is  given  by 

P(S0j|D0j  A  Xq) 

As(Xg)=  P(S0jlDQj  A  Xq) 


(2.10) 


(2.11) 


Note  that  the  denominator  (normalizer)  terms  of  (2.10)  are  relatively  simple  and  are  computed  only 
once  for  each  segment. 
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SECTION  3 


TERMS  OF  EXPECTED  LIKELIHOOD  RATIO 


State  Prior 

P(Xq)  is  the  prior  pdf  of  target  state  (geometric)  with  the  replacements.  As  developed  in  A5, 


P(X^)  =  P(R0)  P(V  P(V7  )  J  (3.1) 

P(Rq)  and  P(0Q)  are  prior  distributions  of  range  and  bearing.  Typically  these  are  assumed 
to  be  uniform  distributions. 


P(V?)  is  the  prior  pdf  of  target  speed  and  heading  where  the  arguments  (V  and  7)  are  computed 

from  the  state  variables  (Ro0QR:0 :)  and  the  inferred  times  (tQtj).  Assuming  a  uniform  heading 
1  1  ‘  1 

distribution  (P(y)  =  — )  and  target  speed  normally  distributed  with  mean  V  and  variance  ay  then, 

2  it 


P(V7) 


(3.2) 


where  V  is  computed  from  X'q  and  the  time  difference  (AtQj)  as  given  in  general  by  the  vector 
equation  (DQ  is  the  displacement  vector  of  sensor  location  at  tQ  relative  to  sensor  location  at  tj) 


V  =  I  D0  ♦  R0  -  Rj  I  /Atoj 


(33) 


The  Jacobian  (J)  transforms  the  Rj0j  variables  into  V7  variables  as  given  by 


_  5l 


V  At; 


(3.4) 


oj 


Thus,  P(Xq)  (3.1)  is  easily  computed  for  the  measured  values  (0Q0j)  and  the  assumed  ranges  (RQRj). 


TIT 


N(.)  is  notation  for  a  normal  distribution  as  given  by 
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f 

To  visualize  the  effect  of  this  prior  state  density  term  (P(Xq))  within  the  RQRj  integral  of  (2.10), 
consider  the  following  as  developed  in  B2.  For  every  V7  value  (given  0o0j)  corresponding  RQRj  values 
can  be  computed, hence  points  in  V7  space  can  be  mapped  into  RQRj  space  as  illustrated  in  Figure  3-1. 

Assume  first  that  V  is  set  at  its  mean  value  (V)  and  7  allowed  to  take  on  all  values  (0  to  2 n ). 

The  resulting  contour  in  RQRj  space  is  an  offset  ellipse  tilted  at  45°  as  illustrated  in  Figure  3.1  (dashed 
ellipse).  The  contours  for  ±0 y  are  also  shown. 

For  the  given  bearings,  the  ellipse  offset  is  proportional  to  own  ship  speed  (VQ)^  and  the  semi 
major /minor  axes  are  proportional  to  target  speed  (V).  Thus  the  la  V7  distribution  (also  shown  in 
Figure  3-1)  maps  into  RQRj  space  as  illustrated.  Within  the  la  region  of  RQRj  space,  P(V7)  will 
be  high,  decreasing  to  negligible  values  elsewhere. 


Figure  3*1.  Target  State  Prior  in  R0Rj  Space 


Now  consider  the  development  of  the  second  term  in  (2.10);  the  detection  term. 

DETECTION 

P(DQj  I  X'q)  is  the  probability  of  detection  occurring  at  times  tQtj.  Neglecting  array  gain  variation 
with  bearing,  then  probability  of  detection  is  a  function  of  range  only,  so  that 

P(Doj  I  X’G)  -  P(D0  I  R0)  PfDj  |  R0Rj)  (3.5a) 

(1)  If  the  (o  j)  measurements  are  made  from  different  platforms  (e.g.,  two  ships  or  offboard  arrays) 
then  offset  is  proportional  to  platform  separation 
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If  more  than  a  single  candidate  is  included  in  the  jth  segment  set,  the  additional  segments 
( t  =  1,2,...)  result  in  an  additional  factor;  the  probability  of  detection  for  these  segments  so  that 

P(Doj  I  XG>  =  p(Do  I  V  P<Dj 1  RoRj>rT  P<D  £  1  R  ,  <XG>  (3.5b) 

where  range  (R  )  is  predicted  from  X'q  at  the  applicable  times  (t  | ). 

As  previously  discussed,  a  target  is  tracked  through  a  convergence  zone  and  smoothed  to  obtain 
an  estimate  of  Yj  (e.g.  bearing,  bearing  rate  and  frequency  rate).  The  time  where  the  putimpfa.  applies 
is  typically  set  at  the  track  midpoint  or  endpoint.  If  this  is  the  only  time  information  used  (total  track  time 
information  ignored)  than  the  P(D|R)  function  would  be  derived  from  the  propagation  loss/signal 
excess/probability  of  detection  curves. 

MEASUREMENTS 

F 

v  is  the  measurement  term  of  (2.10)  as  given  by 


Assuming  the  measurement  errors  are  normally  distributed,  then 

,  .  EL  // 

AY  is  the  difference  between  the  measured  and  true  values;  A  Y^  is  the  transpose  of  AY 
the  conditioning  (Xq)  in  (3.6)  the  true  values  are  predicted  from  Xq. 

W  is  the  inverse  of  the  diagonal  measurement  covariance  matrix  (M„).  M„  has  elements  such  as 
J  y  y 

°eo'  etc'  N°te  that  e°  80(1  ^oue  uncorrelated  with  midpoint  smoothing. 

|Wy|  is  the  determinant  of  Wy 
Ny  is  the  number  of  measurements. 


(3.7) 

.  Due  to 
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As  shown  in  Cl,  the  solution  to  (3.7)  is 

(3.8) 

(3.9) 

AYj  is  the  difference  between  the  measured  values  (excluding  <fo0j)  and  those  predicted  from 
X'G.  Recall  X'G  =  [RoRj0o0j). 

Wq  is  a  diagonal  2x2  inverse  covariance  matrix  with  elements  (1/og^,  1/og^). 

W|  is  a  diagonal  x  inverse  covariance  matrix  with  reciprocal  variance  elements  such  as 
2  .2  2  2 

d/^o*  1^afo*  ^a0j*  l/°fj  ).  Note  the  variables  included  are  all  the  measurables  except  0Q0j. 

Hj  is  the  partial  derivitive  matrix  of  the  measurables  with  respect  to  0Q0j  (de^d^/dd  Qd ^). 

Nj  is  the  number  of  measurables  exceeding  the  two  bearing  measurement,  i.e.,  Nj  =  Ny  (the  total 
number)  -2. 

Note  that  the  calculation  of  W  given  by  (3.9)  is  fairly  simple  since  the  inverse  is  only  2x2. 

Refer  to  C2  for  an  approximation  to  (3.8,  3.9)  and  C3  for  an  example. 

It  should  be  noted  that  the  results  given  by  (3.8)  are  general  and  apply  to  the  association  of  more 
than  a  single  past  segment. 

To  illustrate  the  effect  of  this  Fy  term  within  the  RQRj  integral  of  (2.10),  consider  the  following: 
Assume  that  0  Q  and  fQ  are  measured  in  addition  to  0  Q0y  Using  an  approach  similar  to  that  of  Figure  3-1 
the  la  dQ  fQ  contour  maps  into  RQRj  space  (refer  to  B3)  as  illustrated  in  Figure  3-2. 
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Figiue  3-2.  Rate  Measurements  in  RQRj  Space  —  Track  o 


Figure  3-3  shows  the  la  contour  of  0Qfo  (from  Figure  3-2)  and  also  a  corresponding  contour  for 
the  jth  segment. 
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Figure  3-3.  Rate  Measurements  in  RqRj  Space  —  Tracks  o  and  j 
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The  amount  of  overlap  shown  by  the  shaded  regions  is  related  to  the  magnitude  of  Fy,  the 
conditional  (RQRj)  likelihood  function  of  the  geometric  measurables. 


FREQUENCY  SPECTRUM  LIKELIHOOD  RATIO 


The  frequency  spectrum  likelihood  ratio  (Ag  )  is  developed  in  Appendix  D  and  E  for  two  segments 
with  spectra  SQ,Sj.  In  the  A  formulation  of  (2.10)  the  spectrum  likelihood  ratio  is  conditioned  on  X'G 
and  hence  Ag  should  be  evaluated  for  all  the  RqRj  combinations  in  the  integral  (summation)  of  (2.7). 

The  formulation  in  Appendix  D  assumes  that  the  spectra  have  been  compensated  for  doppler  shift  (a 
function  of  range  rate  difference  (AR)),  where  AR  is  determined  from  Xq.  The  formulation  also  assumes 
that  the  distribution  of  receiver  amplitude  for  any  given  frequency  cell  is  the  same  for  each  segment. 

This  implies  that  the  range  dependence  is  ignored.  A  summary  of  the  results  of  the  frequency  spectrum 
development  is  as  follows. 


The  general  result  is  given  by 


(3.10) 


II  is  product  notation.  n£  is  the  number  of  frequency  cells  that  match.  A  match  is  said  to  occur  when  a 
line  in  segment  o  and  one  in  segment  j  appears  at  the  same  frequency  (after  doppler  correction),  n j  is 
the  number  of  cells  where  there  is  a  mismatch  (a  line  in  only  one  segment).  Ag2  and  Agj  are  functions 


of  the  amplitude  of  the  lines  and  the  X  terms  are  constants  derived  from  the  a  priori  statistics. 


Ag2»  the  “match”  likelihood  ratio,  is  given  by 


(3.11) 
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where 


AS  =  (SQ  -Sj)  =  the  amplitude  difference  (db) 
o  =  amplitude  fluctuation  std  deviation  (dB) 

8S  =  ave  amplitude  relative  to  threshold 

Pp  =  probability  of  fade. 

Agj,  the  “mismatch”  likelihood  ratio,  is  given  by 

SE 

AS1=~  e‘8(j2a)  (3.12) 

where 


SE  =  amplitude  relative  to  the  threshold  (dB). 

These  two  A  terms  are  illustrated  in  Figure  3-4. 

The  X  and  KQ  terms  are  a  function  of  the  a  priori  statistics  of  the  target  and  extraneous  lines  as 
developed  in  Appendix  D,  E.  Definitions  of  these  statistics  are 


nj  =  mean  number  of  target  lines, 
ng  =  mean  number  of  extraneous  lines 
Pp  =  probability  of  fade 

np  =  number  of  cells 

and,  , 

nE  “  nE  +  nT 

Substituting  the  X’s  from  Table  D-l  of  D4  yields 


nO 

(1-Pp)  n-p  n^ 

AS  - 

TTZ 

AS2  (  ,  o  ,+  1 

nE2 

-  . 

AS1  < 


Pp  n»p 


n>j>(l*Pp) 


(3.13) 
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If  amplitude  data  is  not  made  available,  the  spectrum  is  referred  to  as  “clipped”,  then  only  the 

frequency  cells  containing  a  detection  are  available.  For  this  case,  Ag£  and  Agj  are  unity  and  hence 

.  .  n9  *  no 

the  products  (e.g.,  u  (  ) )  can  be  replaced  by  the  exponent  form  (e.g.,  (  )  )  so  that 

!ni 

nT(l-PF) 


(1  -  PF)  nTnc 
-  +  1 

n2 

^FnT+^nE 

n'g2 

2"e' 

Consider  a  typical  case  where  Pp  =  .2,  nq,  =  2,  ng  =  2  and  n^  =  400  (ng'  =  3.8). 
then, 

Ag  =  (44  +  l)"2  (0.05  +  0.53)ni  6 


(3.14) 


or 


(3.15) 


Log  Ag  =  1.65  n£  -  0-24  nj  +  0.78 


This  function  is  plotted  in  Figure  3-5. 

Another  approach  was  initially  used  to  develop  the  spectrum  likelihood  ratio  (Appendix  F).  It  was 
dropped  in  favor  of  the  approach  of  Appendix  D  basically  for  two  reasons;  amplitude  information  is  not 
utilized  (clipping  is  assumed)  and  no  extraneous  lines  are  included. 


DENOMINATOR  TERMS  (Reference  (2.10)) 

•  P(R0  )  —  State  pdf  R  and  6  are  considered  essentially  independent  so  that 


P(Ri®j[)  *  P(Rj)  P«V 


(3.16) 


•  P(D|R)—  Same  as  numerator  detection  terms. 

•  P(flflDR)  —  As  shown  in  A4,  this  pdf  is  obtained  by  including  relative  cross  range  rate  (X)  and 

integrating  it  out  as  given  by 


P(0f|DR)  =  /p(ef|DRX)  P(X)  dX 
X 
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n1  *  NO.  MISMATCHES 


Figure  3-5.  Log  Likelihood  Ratio  for  Frequency  Spectra 


..  2  2 

Assuming  normally  distributed  8  f  errors  with  variances  (o'q  oj  )  then, 

,  r  v2 . 


(o  -X/R\  Af-X2/R\ 


(3.17) 


Assuming  normally  distributed  a  priori  target  speed  with  mean  V  and  variance  oy  and  uniformly 
distributed  target  heading  then, 


Iff  V  -V,  -  -1 . -~1 

P(X)~  it  v=|Xt|LN(  oy  JV2-Xt2 


dV 


(3.18) 


where  Xj  is  target  cross  range  rate  (Xj  =  X  +  XQ),  XQ  is  own-ship  cross  range  rate. 

This  function  is  shown  in  Figure  3-6  for  V  =  15  kts  for  XQ  =  0.  For  finite  XQ,  simply  shift  the 
curves  left  by  an  amount  XQ. 


Figure  3-6.  Cross  Range  Rate  Prior  Distribution 
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SECTION  4 


LIKELIHOOD  SUMMARY 


This  basic  A  expression  and  the  terms  included  are  summarized  below.  (Refer  to  Appendix  A.) 
//  P(Xq)  •  P(D0j  I  X'G)  •  Fy  •  AS(X'G)  dR0Rj 


_  R0Rj 


n  ^  PlR^j)  •  P(Dj  I  Rj )  •  P(Yl4  I  DjRj)  dRj 
Target  State  Prior 

p(x'G)  =  p(RQ) '  p(«o)  *  pv7J 


(2.10) 


(3.1) 


Detection 

p<Doj|XG)  =  P(DoDj  I  RoRj>  ‘  np(°.  I  XG>  (3.5b) 

Measurements  (Geometry) 

Fy  =  Kyi  e  AYlTwAYl  (3.8) 

Measurement  (Spectra) 

,  n-i  nrp 

AS(XG)=TT  IaS2  X21  +  1lTT1  l^Sl  X11  +  X12l  *  e  1  *  (3.10) 

Mapping  of  the  a  priori  speed/heading,  detection  and  measurement  densities  into  range  (RQRj) 
space  as  previously  discussed  is  shown  as  a  composite  in  Figure  4-1.  Note  that  bearing  rate  and 
frequency  rate  are  the  only  measurements  included  here.  Normalization  (A^),  the  Rj  factor  of 
(3.4)  and  the  frequency  spectrum  term  are  not  included  in  this  illustration. 

The  association  measure  A  is  proportional  to  the  summation,  over  the  RQRj  space,  of  the  product 
of  terms  in  (2.10).  Thus  only  the  regions  of  overlap  of  all  three  terms  of  Figure  4-1  will  contribute  to 
the  sum.  Note  that  the  primary  contribution  is  from  the  zones  at  RQR-  3  60,90;  the  true  range  pair. 


A  block  diagram  showing  the  process  used  to  compute  AQj  and  F0j  is  shown  in  Figure  4-2. 
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MEASUREMENT 
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Figure  4-1.  Composite  Mapping  into  RQRj  Space 
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-  ASSOCIATION  MEASURE  AND  LIKELIHOOD  RATIO  DEVELOPMEN’ 


A1  -  DERIVATION  OF  ASSOCIATION  MEASURE  P(Aoj  |  Zojk) 


The  general  problem  is  to  determine  how  well  a  current/reference  segment  (subscript  o)  associates 
with  one  or  more  candidate  segments  (subscript  j).  The  remaining  segments  are  identified  with  sub¬ 
script  k  so  that  the  total  data  set  is  defined  by  ZQjk.  Thus  the  probability  of  association  of  the  current 
(o)  and  candidates  (j)  given  all  the  data  is  expressed  by 

Foj  =  p(Aoj  I  20jk)  (1) 

Note  that  the  j  and  k  segment  sets  consist  of  one  or  more  segments. 

Applying  Bayes’  rule  yields 


p,.  IZ  , 

'  °j 1  Z°jk  '  P(Zojk) 


Consider  the  first  numerator  term.  The  k  subscripted  segments  are  not  conditioned  by  Ac 


P^ojk  *  A0j)  P(Zoj  I  A0j)  • 


Consider  the  following  simple  equality 


P(ZQ)  •  [P(Zn)  •  P(Zj2) ... ) 

P,Zk) '  • p,z*’ 


Note  that  the  numerator  is 


P(ZQ«  Zj,  Zj^  I  A0jk)  A  P(Z0jk  I  AQjk) 


Where  A0jk  is  the  condition  that  the  ojk  (i.e.,  all)  segments  are  not  associated. 


Also  note  that  the  denominator  is 


For  m  past  contacts  of  which  none  or  one  are  associated  with  the  current  (o)  segment  then 


m 

P(Z)  =  P(Z  |  A)  P(A)  +  2  P(Z  |  Aoi)  P(Aoi) 
i 


Substitution  into  (11)  yields 


p(Aoj  I  Z)  = 


oj 


P(A)  +  ”  P(Z  I  Aoi)  P(Aoi) 


P(A°j)  >  P(Z|A)P(Aoj) 


Note  that 


(14) 


(15) 


P(Z  I  Aoi) 


p<ZoilAoi> 


P(Z|  A) 
P(Zoi'Aoi) 


=  AoiP(Z|  A) 


(16) 


Assume  that  all  P(AQi)  are  equal  and  substitute  (16)  into  (15)  to  obtain 


p(Aoj|Z)  = 


oj 


P(A) 

p(Aoj) 


m 


+  2  ^i 


(17) 


Assuming  a  diffuse  prior  where  Pg  is  the  probability  that  the  correctly  associated  past  segment 
exists  in  the  set  of  m  segments.  Then  P(A)  is  the  probability  that  it  does  not  exist  as  given  by  (1-Pg). 
P  (A0j)  is  the  probability  that  the  j  th  segment  is  associated  with  the  o  segment  as  given  by  Pg  •  1/m 
so  that  the  first  denominator  term  is  ((1-Pg)  /  pg)  m.  Thus, 


p(Aoj  I  Z)  = 


(18) 


A2  -  A  FORMULATION 


The  oj  subscript  on  A  will  be  dropped  for  simplicity,  from  Al(12) 

_  P(Zoj'Aoj) 

P<Zoj  I  V 


(1) 


A -3 


Since  segments  o  and  j  are  unassociated  (denominator)  they  are  independent  so  that 


The  measurement  vector  consists  of  detection  (D),  geometric  (Y)  and  spectrum  (S)  elements 
(ZQj  =  (DYS]0j)  (refer  to  Section  2  (1.1)).  Expanding  the  first  intergrand  term  yields 

P(ZQj|X)=  P(DYS|X)  (2) 

=  P(D  |  XGXg)  P(Y  I  DXgXs)  P(S  |  DXgXs) 

where  subscripts  (oj)  apply  for  DY  and  S.  The  detection  event  is  conditioned  on  source  spectra  (Xg), 
a  fairly  complex  relationship  hence  the  Xg  conditioning  will  be  removed  from  the  first  term.  The  Xg 
condition  is  irrelevent  in  the  measurement  (2nd)  term  and  hence  will  be  dropped. 


A4 


Then,  substitution  into  (1)  yields 


The  geometric  state  variables  selected  are  range  and  bearing  for  both  the  current  (o)  and  the  most  recent 
past  segment  (j)  so  that 


XG  =  [VoW  (4) 

where  the  corresponding  times  (tQtj)  are  implied. 

Bearings  (0Q0j)  are  always  measured.  (~)  will  be  used  to  identify  measured  values. 

Thus,  from  (3) 

An  =  /  P(XG)  P(D|XG)  P(Y|DXG)  P(S|DXG)  dXG  (5) 

XG 

Since  P(XG)  is  a  weak  function  of  0Q0j,  a  good  approximation^  is  to  assume  that  the  state 
variables  0Q0j  equal  the  measured  values  (djs A  Prime  (')  will  be  used  to  indicate  this  substitution, 
i.e.,  P(XG)  ♦  P(XG).  The  same  applies  to  the  detection  and  spectrum  density  functions.  Thus  all  except 
the  measurement  density  P(Y|DXG)  can  be  moved  outside  the  0Q0j  integral  so  that 


<  )  This  is  a  good  approximation  for  temporal  association  where  time  between  segments  is  large 

(>  30  minutes)  as  it  is  for  segments  in  different  convergence  zones.  For  short  times  another 
formulation  applies  (Refer  to  A6). 
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A4  -  A  d  DEVELOPMENT 


A,j  is  defined  as  the  denominator  of  A  (Refer  to  (2.4a).  The  state  variables  selected  are  range, 
bearing,  range  rate  and  cross  range  rate.  Thus  for  a  typical  o  th  segment 

XG  =  Ro0oRoXo  (1) 

The  spectrum  component  of  the  state  vector  is  Xg,  thus  for  a  typical  segment  (similar  to  A3(3)) 

Ad0  =  //  P(D0IXG)P(Y0ID0XG)P(S0ID0XGXg)P(XG)P(Xs)dXGXs  (2) 

XGXS 

Assuming  detection  is  a  weak  function  of  bearing  (omni  array  gain)  then 


P(D0IxG>  =  P(D0'R)  <3> 

~  ~  -v 

Assuming  that  the  measurables  are  bearing,  bearing  rate  and  frequency  rate  (0Q0  0fQ)then,  with 
o  subscripts  implied, 

P(Y  |  D  XQ)  *  P(0  I  D  0 )  P{0  f  |  D  RX)  (4) 

Note  that  R  is  not  included  since  0  f  are  a  function  of  R  and  X  only;  not  R. 

The  spectrum  is  only  a  function  of  range;  doppler  shift  for  a  single  segment  is  irrelevent  so  that 
P(S  I  D  XG)  =  P(S  I  D  R)  (5) 

The  state  variables  are  essentially  independent  so  that 

P(XG)  =  P(R)P(0)P(R)P(X)  (6) 
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Substitution  for  these  four  functions  into  (2)  and  separation  of  integrals  yields 


Ad  =  [/P(»lfl)P(fl)d»] 

°  e 

j  /  P(D  |  R)  [  /  P(S  |  D  RXg)  P(XS)  dXg]  (7) 

| R  XS 

•l/P(d  fl  DRX)P(X)dXl  P(R)dR[  •  / P(R)dR 
X  )  P- 

* - - -  P(0flDR  — J  L_unity-J 

The  first  integral  reduces  to  the  value  of  P(0Q)  (bearing  prior)  evaluated  at  the  measured  value 
(0=0)  and  the  last  is  unity.  Thus  with  the  indicated  definitions  and  ignoring  the  range  dependence 
(condition)  of  the  spectrum  density 

i 

A d0  *  p(*~o>P(So'Do>  /P(D0'Ro>P<W  D0R0>P(Ro>dRo  |  (8) 

° _  Ro _ I 

The  expression  for  other  segments  has  the  same  form,  simply  replace  subscript  o  with  the  appropriate 
segment  subscript. 

A5  -  PRIOR  STATE  -  P(Xq) 

The  prior  pdf  of  the  geometric  range,  bearing  state  vector  at  the  time  of  segments  o  and  j  is 
given  by 


p<xb>  =  P(VoRjV 


(1) 


Range  and  bearing  densities  are  essentially  independent  so  that 

p(x'G)  =  p(rd)  p<r0)  p(Rjrj  i  R0ro) 


(2) 
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The  third  term  can  be  evaluated  from  the  known  target  speed  and  heading  density  P  (V  7)  using  the 
Jacobian  (J)  transformation  as  given  by 


P(Rjflj  |  R0(T o)  =  P(V  7)  ’  J  (3) 

where 


av7 

>1 

<©  1 

>1 
r©  1 

_ 

cos  0  R; 

J  sin  0 

3Rj0j 

3Rj  30 j 

3t_  3t_ 
3Rj  30 j 

At 

sin  0  Rj 

COS  0 

VAt  VAt 

where 

I  |  is  determinant  notation. 

0  is  target  heading  relative  to  the  line  of  bearing  (0j). 

thus, 


Rj 


V  At^ 


Note  that  target  speed  (V)  is  a  function  of  )  as  given  by 


=  'Do  +  Ro  “  Rj  I  / 


(5) 


(6) 


where 


II  is  magnitude  notation  and  R  is  vector  notation.  DQ  is  the  displacement  vector  of  sensor 

location  at  t  relative  to  sensor  location  at  t-  At  =  t  -  t. 

J  o  y 
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A6-  SHORT  TERM  ASSOCIATION 


In  the  previous  analysis  the  state  prior  could  be  removed  from  inside  the  0O  0j  integral  since  the 
time  difference  between  segments  was  assumed  to  be  large.  This  was  an  important  simplification 
since  the  9  integrand  was  then  gaussian  and  hence  readily  integrable.  This  assumption  is  valid  for 
CZ  to  CZ  temporal  association  however  it  breaks  down  for  short  term  or  for  mutual  (simultaneous) 
association. 

This  analysis  develops  the  short  term  likelihood  expressions  using  a  different  state  vector. 

In  the  long  term  temporal  association  formulation  Xq  was  defined  as  the  range/bearing  at  the 
times  of  the  two  segments  as  given  by  A3  (4). 

For  short  term  association  a  more  natural  state  vector  is 


XG  -  l\eoVy] 

"K  O 

where  V7  are  target  speed  and  heading;  V  is  assumed  to  be  gaussian  with  mean  V  and  variance  Oy  . 
Then  A3(5)  can  be  written  as 


//  P(D|Ro^^7)P(S|DX^)P(^6T7)[//P(V)P(Y|XG)d0oV]dRo7  (1) 


V 


X'A 


F' 

*  IF 


Note  that  the  detection  and  spectrum  densities  have  been  removed  from  the  9  QV  integral  since  they 
are  weak  functions  of  O^and  V  and  hence  the  measured  value  of  90(9Q)  and  the  a  priori  (mean)  value 
of  target  velocity  (V  )  are  used. 
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For  comparison  the  previous  0QOj  integral  and  the  0QV  integral  are  rewritten  below 

•  Longterm  (refer  to  Cl(l)) 

Fv  -  f/P(YlXG)d0o0rffP(Ool0o)- (2) 

°j  Vj 

•  Short  term 

Fy'  =  P(V|V)P(Y|XG)d0oV  =  ^P(e'o|0o)P(V|V)P(YilXG)d0oV  (3) 

0oV  *oV 

Note  that  the  forms  are  the  same  except  that  0j  is  replaced  by  V  and  that  Yj'  includes  all 
measurables  except  0  Q  (recall  Yj  excluded  both  8  Q  and  0j). 


Referring  to  (3.6)  and  its  solution  (3.8)  the  solution  to  the  0QV  integral  is  given  by 


i 


IW'I 


(2ir)Nl 


e->MYi  W’AYj 


(4) 


where 


W'  =  Wj  -(W'jH'j)  (WJ,  +  HiT  W'j  H')1  (W'j  H'j) 


(5) 


AYj'  is  the  difference  between  the  measured  values  (excluding  0Q)  and  those  predicted  from 


XG" 


where 


XG  ’  “  IVoV>l 


W^'  is  a  diagonal  2x2  inverse  covariance  matrix  with  elements  (1  log^,  l/o^). 

Wj'  is  a  diagonal  Nj'  x  Nj'  inverse  covariance  matrix  with  elements  such  as  (l/o^,?  l/o'g 
Note  the  variables  included  are  all  the  measurables  except  0Q. 

H  {  is  the  partial  derivitive  matrix  of  the  Yj  measurables  with  respect  to  0  Q  and  V . 


Nj  is  the  number  of  measurables  excluding  0Q.  Thus  if  Ny  is  the  total  number  Nj  =  Ny-1. 
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For  comparison,  the  previous  (long  term)  numerator  likelihood  ratios  and  the  short  term 
expression  (1)  are  rewritten  below 

•  Long  Term 

An  =  / /  P(Wjf)-  P(D|Xg')P(S|DXg')-  Fy  dR^  (6) 

R0Rj  - 

XG 

•  Short  Term 

An  =  ffp(R o*0y) '  P(dIXg")  *  P(S|DXg”)  •  F  '  dRQ7  (7) 

V 

The  variables  of  the  prior  term  are  essentially  independent  so  that 

p(RoV)=  p(n0>p<Oph>  <«> 

Again  assuming  probability  of  detection  to  be  a  soft  function  of  azimuth  then 

P(D|Xg'>  P(Dq|R0)  P(Dj|Rj(XG''))  (9) 

The  spectrum  term  is  a  function  of  the  change  in  range  rate  between  the  two  segments  so  that 

P(S|DXg")  =  P(S|DAR(Xg''))  (10) 

Thus  the  short  term  expression  is 

An  =  P(0O)  /P(D0IR0)P(R0)/P(Dj|XG'')P(T)P(S|DXG")[Fy’JdR0r  j  (U) 

R0  7  j 

The  denominator  (A^ )  is  the  same  fo.  the  short  term  and  long  term  formulations. 
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APPENDIX  B  -  GEOMETRY  AND  RQRj  CONTOURS 

Bl-flf  MODEL 

Nonlinear  Own  Ship  Motion /Two  Platforms 

Consider  the  general  case  where  the  o  and  j  th  measurements  are  made  from  different  platforms 
(2  ships)  or  from  the  same  ship  with  displacement  caused  by  own  ship  motion.  In  either  case  DQ  is 
the  displacement  between  the  position  of  the  two  platforms  (ships)  for  the  time  interval  At  *  (tQ  -  tj) 
and  0  is  the  “heading”  of  the  displacement  line  (vector). 

The  geometry  is  shown  in  Figure  B1 . 


81-587 


V  At 


TARGET 


Figure  B1  .  Geometry 
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Definition  of  terms: 


^  t  ~  lj  *  time  interval  between  measurements 
o’  own  ship  and  target  displacements  during  At 
<t>  =  “heading”  of  own  ship  displacement  line/vector 

7,  V  -  target  course  (heading)  and  speed 
Voo’  ro=  own  ship  speed  and  heading  at  time  tQ 
Voj.  rj  =  own  ship  speed  and  heading  at  time  tj 
0 o’  °i  =  true  bearing  at  times  tQ  and  tj 
Ad  =  gQ-  0.  =  bearing  difference 
R0-  Rj  =  Range  at  tQ,  tj 


Solving  for  the  relative  bearing  rates  at  tQ  and  tj  yields 
Rj  sin  Ad  -  DQ  sin  (dQ  -  0) 


R. 


R; 


At 


+  v00  sin  (0O  -  r0) 


(i) 


RQ  sin  Ad  -  DQ  sin  (0 j  -  0) 

^  “  ♦  VQj  sin  «?j  -  Tj) 


Range  acceleration  (R)  equals  (_Af )  where  A  is  wavelength.  R  also  equals  d* 
is  given  by 


R  so  that  frequency  rate 


At 


and 


Xti 


(2) 
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Linear  Own  Ship  Motion 


Assuming  own  ship  is  on  a  constant  course  and  speed  between  points  o  and  j,  then  £  *  rQ  *  Tj 
and  DQ/At  *  VQO  -  VQj.  Then  (1)  and  (2)  reduce  to 


lj  sin  Ad 


At 


«  V 

*  ^  =  r 


o  sin  Ad 


At 


Xt 


R: 


sin  Ad 
At 


R, 


V  -TT 


(3) 


sin  Ad 
At 


B2  -  TARGET  SPEED  CONTOURS  ON  A  RQRj  PLOT 

The  problem  here  is  to  determine  how  the  prior  target  speed/heading  uncertainty  area  maps  into 
R0Rj  space.  The  inputs/givens  are  own  ship  displacement  (DQ  <t>),  the  two  bearings  (dQdj)  and  At. 

Thus  all  lines  in  Figure  B1  are  known /fixed  except  the  target  displacement  vector  (D-p-y)  or  (VAt.y). 
For  a  fixed /assumed  value  of  target  speed  (D-p/At)  variations  in  target  heading  will  generate  the  desired 
range  pairs  (RQRj). 

The  approach  is  to  solve  the  two  triangles  of  Figure  B2  (also  refer  to  Figure  B1 ). 

The  range  pair  at  the  bearing  crossover  point  (T)  is  obtained  using  the  law  of  sines  so  that 


Rso  =  D0sin  (dj  -  0)/sin  Ad 
R8j  -  Dosin(do-0)/sinAd 


(1) 


Note  in  passing  that  this  would  be  the  crossfix  range  solution(s)  assuming  zero  target  speed. 


The  true  range(s)  (RQRj)  differ  from  this  crossover  range  (R^R^)  by  the  increment  ARq  ARj 
due  to  target  motion  as  shown.  Using  the  law  of  cosines  yields  the  relationship  between  ARQ  and 
ARj  given  the  Ad  measurement  and  the  assumed  target  speed  /displacement  (D<p)  so  that 


*  ARq2  ♦  ARj2  -  2  ARq  ARj  cos  Ad 


(2) 
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Thus  given  an  assumed  target  speed  and  At  then  D«p  is  known  and  A  Rj  can  be  solved  as  a  function 
of  ARQ.  The  result  is  a  contour  in  R0Rj  space  as  illustrated  in  Figure  B3.  Rotating  45°  to  a  new 
coordinate  system  AR0'  AR^'  as  shown  and  using  (2)  yields 


1 


ARj 


-2 


V 


(3) 


Note  that  this  is  the  equation  of  a  ellipse  with  semi  major  and  semi  minor  axes  as  given  by 


a0  ~  DT  /  VI-  cos  A0 

(4) 

aj  =  ftj.  /  V  1  +  cos  A0 

Expressing  the  offset  point  (  (l)  )  in  this  rotated  frame  yields 
,  Do 

Rso  =  ~ -  (sin  eQ  +  sin  0:) 

V2  sin  AS  J 

R  '  °° 

Rsj  =  - -  (sm«o 

v2  sin  Ad  J 


(5) 


B3-0  f  CONTOURS  ON  A  RQRj  PLOT 


8  Contours 


Rearrangement  of  Bl(3)  yields  the  RQRj  contours  for  8Q  as  given  by 


R; 


Similarity  for  8 


6QAt 


sin  A6 


Oj  At 
sin  A  8 


Thus  the  8  q8j  contours  are  simply  straight  lines  as  shown  in  Figure  B4(a). 


(1) 


(2) 


f  Contours 


Rearrangement  of  B1  (3)  yields 


Similarly  for  f  j 


(3) 


(4) 


These  f  contours  (3)  and  (4)  are  illustrated  in  Figure  B4(b)  including  a  perturbation  on  fQ  and  fj 
to  illustrate  the  effect  of  f  changes/errors. 


B4  TYPICAL  GEOMETRY 

Consider  the  following  conditions:  RQ  =  60  nmi,  Rj  =  90  nmi,  8 Q  =  45°,  VQ  *  20  kts,  and 
V  =15  kts.  Headings  are  0°  and  180°  for  own  ship  and  target  respectively 
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Using  the  law  of  sines  [60/sin  0  =  90/sin  135°  ]  then  <t>  ~  28.1  ,  thus  AO  -  16.9' 

Using  the  law  of  cosines,  relative  motion  is  obtained  as  given  by 

D  =  [602  +  902  -  2-  60-  90  cos  16.9°  1 1/4  =  37  nmi 
then  the  time  difference  is  given  by  At  =  D/(  V0  +  V )  =  1.03  hr. 


"o 

8  CONTOURS 


R« 


f  CONTOURS 
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Figure  B4.  0  f  Contours 


1 

1 

.  APPENDIX  C  -  MEASUREMENT  LIKELIHOOD  FUNCTION  DEVELOPMENT  -  (Fy) 

1 

Cl- Fy  FORMULATION 

The  measurement  likelihood  function  with  0 Q0  j  integrated  out  is  (refer  to  A3(6)) 

Fy  ■4{P(YoiIDoiXG)d».»i 

(1) 

Assuming  normally  distributed  measurement  errors  and  0Q0j  replaced  by  AXjfor  generality)  then 

F  =  K  J  e'^dAX0 

AX0 

(2) 

where  |Wy|  is  the  determinant  of  the  inverse  of  the  measurement  covariance  matrix  (diagonal)  and 

Q  =  (AYq-HAX0)T  (AYq-HAX0) 

(3a) 

i 

\ 

H  =  3Y/3X 

(3b) 

Ny  =  number  of  measurements  (elements  in  Y  vector) 

(3c) 

i 

and, 

AYo  =Y-Y0(X0),  AXq  =  X-X0 

(3d) 

! 

Ky  =  (2»r)-Ny/2(\nW^I) 

(3e) 

! 

) 

Note  that  the  measurements  (Y)  and  the  partials  (H)  have  been  normalized  with  respect  to  their 

standard  deviations  (o ).  Define  AX'  to  be  the  change  in  X  relative  to  the  optimum  but  unknown 

solution  (X*),  Thus, 

1 

1 

AXq  =  AXq*  +  AX’ 

(4) 

1 

1 

The  optimum/mode  is  given  by 


AXq*  =  [(H1^)'1  HT]  AYq  (5) 

< -  X - ' 

The  transformation  matrix  [  ]  is  designated  X  as  indicated.  Substituting  (4)  and  (5)  into  (3a) 

and  expanding  yields 

Q  =  AY0T[(I-HX)T(I-HX)]  AYq  (6) 

-2aX'T[HT(I-HX)]  AYq 
+  AX'T  [  HT  H  ]  AX' 

Note 

(1)  [  (I -HX)T  (I  -  HX)  ]  reduces  to  [I-  HX]  after  noting  that  XH  is  an  identity  matrix 
and  HX  is  symetric. 

(2)  [Ht  (I  -  HX)  ]  is  zero  since  expansion  yields  [HT  -(H*H)  (HTH)‘*  HT]  and 

(  )(  r1-!. 


where  Kyl  =  Ky/k  and  Nxis  the  number  of  variables  in  X'  (number  integrated  out). 

Using  the  definitions  of  Qy  and  X  in  (7)  and  (5)  yields 

Qy  =  A  Yqt  [  I  -  H  (H^H)'1  HT  J  A  Yq  (10) 

Recall  that  A  YQ  and  H  are  measurements  and  partials  that  have  been  normalized  with  respect  to 
oy.  Define  the  inverse  covariance  matrix  of  Y  to  be  Wy.  This  is  a  diagonal  matrix  with  elements  ay'2. 
Define  the  matrix  >/  Wy  to  be  the  matrix  obtained  by  taking  the  V  of  each  element.  Using  the  (A) 
subscript  to  indicate  unnormalized  (actual)  measurements  and/or  partials  yields 


Partition  A  Y^  into  the  variables  integrated  out  (YQ)  and  those  remaining  (Yj).  WQ  and  Wj  are  the 
corresponding  inverse  covariance  matrices.  Thus 


Substitution  into  (9)  and  (11)  and  matrix  manipulation  yields  Qy  and  W  as  given  by 

Qy  =  AY1T(W1-(W1H1)(W0+H1TW1H1r1(W1H1)T]  AYj  (13) 

I -  W  - > 

Ky  Factor 

The  Ky  factor  (9)  could  be  determined  by  matrix  algebra,  however  a  more  straightforward  approach 
is  as  follows.  Note  that  the  density  function  integrand  of  (1)  is  assumed  normal  so  that  Fy  is  also  normal. 
From  (9)  and  (13) 
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(14) 


Fy  =  Kyle-**Y1  ' ^Y1 


In  general,  the  n  dimentional  normal  distribution  of  x  with  inverse  variance  Wx  is  given  by 


|WX|  e‘^xTwxx 


By  comparison  of  these  two  normal  distributions  then, 


(2tt)N1 


where  Nj  is  the  number  of  measurables  in  Yj  as  given  by  (Ny  -  Nx). 


Final  Solution 


Substitution  of  (16)  and  (13)  into  (9)  then  yields 


IWj -  .-KAYiTwAYi 


(2ir  )N1 


where. 


W  =  Wj  -(WjHj)  (WQ  +  H^WjHj)-1  (WiMj)1, 


C2-  APPROXIMATE  SOLUTION 


Consider  the  use  of  l/o  subscript  notation  to  indicate  effect/cause  respectively.  Thus  Wj  as 
currently  written  would  be  replaced  by  (W  with  measurement  of  type  1  was  caused  by 
measurement  type  1  (same)).  Wj^0  would  then  mean  W  with  measurements  of  type  1  caused  by 
measurements  of  type  o.  Thus  [  ]  of  (13)  is  rewritten  as 


w-ll-  Wi/i-Wi/i[Hi/0(Hi/0TW1/1H1/0^W0/0)-1h1/0Tjw1/1 


A  covariance  matrix  will  be  designated  by  M,  the  inverse  of  W  (i.e.,  It  can  be  shown 

that  a  covariance  matrix  of  one  variable  type  can  be  transformed  to  another  by  the  matrix  operation 
HMHT.  For  example,  transforming  from  type  o  to  type  1,  yields 

Ml/o  “  Hl/o  ^o/o  Hi/o  (2 

T 

Similarly  the  inverse  covariance  matrices  can  be  transformed  using  the  operation  H  WH.  Thus 


Wo/l  =  Hl/o  Wl/1  Hl/o 


then  the  [  ]  j  term  of  (19)  can  be  dissected  as  follows 


■  Hl,o  <H1/oT  »1,1  H1/o  *  W1  H1/0T 


*o/o,l  - 

—  ^0/0,1 


Thus  [  ]  j  =  i  is  the  joint  covariance  matrix  of  type  1  variables  after  estimation/smoothing  using 
both  type  1  and  o  variables.  It  is  obtained  by  determining  the  smoothed  covariance  matrix  in  terms  of 
type  o  variables  (M0^0  j)  and  then  transforming  to  type  1  variables. 

An  alternate  approach  to  obtaining  j  and  easier  to  interpret  and  simplify,  is  to  first  transform 
W 0j0  to  Wj^0  and  then  combine/smooth  in  type  1  variables  so  that 


-  Ml/o,l  •  I  <Hl/o  “W1  H,,/)-1  *  wl,l  l'1 


u^o/oJ 
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Note  that  this  would  simplify  if  M^Qwas  assumed  to  be  diagonal,  then  W^0  and  [  J^are  diagonal 
and  (19)  reduces  to 

W*  =  Wl/1  -  W1A  (W1/G  ♦  W^)1  W1A  (24) 


Since  is  diagonal  and  W^Q  assumed  diagonal  then  W*  is  also  diagonal.  Thus  any/each  diagonal 
element  can  be  determined  independent  of  the  others  and  inverses  become  reciprocals.  Consider  the  i  th 
diagonal  element  and  for  simplicity  let  a  =  Wjyj  and  b  =  ^\/0  . 


Then, 


W: 


ai-aj 


b-  +  fl* 
ui  1 


(25) 


which  reduces  to 
Wi* 


1/a*  +  1/bj 


(26) 


o 

where  it  is  noted  that  1/aj  and  1/bj  are  the  covariances  and  Mjy0j  which  are  simpb'  and 

°l/oi^-  Note  is  the  variance  of  a  type  1  measurement  (Yj).  a  i/G2  is  the  variance  of  the  errors 

in  predicting  Yj  due  to  the  type  o  errors  (YQ). 


Mi/o  j  is  obtained  using  the  transformation  of  (20)  so  that 


Wj*  =  - - ~  (27) 

‘^1/1  i  +  Hl/o  i  ^o/o  ^l/o  i 


For  notational  convenience  this  will  be  written  as 


W.*  =  - - - —  (28) 

Mii  +  Hli  M0  HjT 

Recall  that  Mjj  is  the  i  th  covariance  (ogQ^,  ogf,  ...),  (6Q6 j  not  included).  Hjj  is  the(Nj  x  2)  partial 

derivitive  matrix  of  the  Nj  type  1  variables  with  respect  to  the  two  type  o  variables  (0o0j).  MQ  is  a 

2  2 

2x2  diagonal  covariance  matrix  with  diagonal  elements  OgQ  ■ 


CS 


Thus,  since  W*  is  diagonal,  replacing  W  by  W*  in  (13)  yields  the  simple  expression 


~  — 

-  — 

-  r0' 

H11 

1'. 

-  r  0  ,  - 

H12 

+  2r  Rj00'  ,  - 

A 

H13 

+  2r1Ro00',  - 

H14 

-  _ 

. _  _ 

9  <Wof  j 
3flo0j 


where  X  is  wavelength,  Xf  is  considered  the  measurable  rather  than  f  and 


r 


A  •  A  sin  Afl 
Ro’0  At 


A  cos  AO 
At 


(2) 


(3) 


AO  and  At  are  the  measured  bearing  and  time  differences.  The  second  column  of  Hj  is  the 
negative  of  the  first  column. 


Wj  is  4  x  4  (Nj  -  4),  diagonal  with  elements  that  are  the  reciprocal  of  the  measurement  covariances 
(Mj)  elements  0XfQ^  °\f2^  (  )  are  the  elements  of  Mjj,  i  =  1,2,3, 4.  WQ  (similarly)  is  2  x  2 

diagonal,  with  elements  that  are  the  reciprocal  of  the  0 Q0  j  covariances  (MQ)  elements.  If  the  simplified 
form  of  Wj*  is  used  then  the  Mj*  covariances  are 


i 

Mj*  =  1/Wj* 

1 

<’»o2  +  N')2»02 

2 

ag.2  +  (r1  0')2  0o2 

3 

oXfo2  +  ( 2rRj00')2  a  2 

4 

axfj2  +  (2r-1Ro00,)2ao2 

where, 


7o2~  °0O2  +  °0j2 


then,  for  example,  one  of  the  Fy  factors(F y^Hepplying  to  the  0  ^  measurement)  is  given  by 

1 


yi 


exp 


[-p- &->'*,■] 


(4) 
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=  V  2n  Mx* 


(5) 


C4  -  REDUCTION  USING  A6  FOR  SINGLE  SHIP  WITH  LINEAR  MOTION 


As  shown  in  Bl(3)  the  measurables  0ofo®jfj  (plus  8Qd-)  are  a  function  of  Ad  (and  RqRj)  and  hence 
the  double  integral  of  Cl  (1)  reduces  to  a  single  integral  with  variate  AO  having  variance  aA0  2  = 

°$  +  °6^4  Thus  the  term  (  )'^  of  Cl(13)  is  a  scalar  given  by  (note  Hj  =  3Y^j/9A0) 


Thus  manipulation  of  Cl(13)  yields 


(1) 


w  -  W1-(W1H1)(W1H1)TMo/o  I 
Q  ■aV1Tw1iY,-(H1TW14V1,2Mo/ol 


/*Yli\2 

v  u 

V*u  /  " 

1  2 
°Yli 

(2) 

(3) 
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APPENDIX  D  -  DEVELOPMENT  OF  FREQUENCY  SPECTRUM  LIKELIHOOD  RATIO 


Dl~  LIKELIHOOD  RATIO 


Assuming  cell  to  cell  independence  and  n  cells  then 


"c 

=  "  A;  (1) 


where  the  likelihood  ratio  for  a  typical  ith  cell  is  given  by 
P(S:Dj|A) 

Aj  =  - -  (2) 

PiSjDjIA) 

SjD-  is  the  amplitude/detection  information  for  the  ith  frequency  cell.  A  is  the  condition  that  the  two 
data  segments  (complete  frequency  spectra)  are  associated;  i.e.,  both  have  a  common  origin  (target). 

A  means  they  are  not  associated.  These  can  be  interpreted  as  0-1  hypothesis  (Hj  -  A,  HQ  =  A). 

D2— MEASUREMENT  AND  STATES  DEFINED 

The  target  will  be  modeled  as  a  source  with  n^  frequency  lines.  This  is  the  mean  number  of  lines 

detected  in  two  segments  (segment  oj).  The  line  amplitude  (Xj)  (without  fluctuation)  is  defined  by  the 

probability  density  function  P(Xj)  or  simply  P(X).  Fluctuation  is  considered  log  normal  with  standard 
2 

deviation  o^.  Probability  of  fade  is  Pp.  Extraneous  lines  are  also  present;  uncorrelated  from  segment 
to  segment  with  the  mean  number  given  by  ng. 


Thus  the  measurements  are 


r—  n 

So 

”soi- 

.  soi.. 

•  Son 

_si_ 

_sji  ■ 

.  Sji .. 

■  v  J 

-  segment  o 

-  segment  j 


and  the  target  variables  are  (nT,  op,  Pp)  and  amplitude  Sx  with  Pdf  P(SX). 


The  prior  conditions  are  defined  in  terms  of  detection  (D)  and  correlation  (C)  states  for  each 
association  condition  (A  or  A). 


D-l 


There  are  four  binary  detection  states  Djj,  Dqj,  Djq,  Dqq  where  Djj  means  the  target  (T)  is 
detected  in  both  segments,  Dqj  means  no  detection  in  one  segment  and  detection  in  the  other,  etc. 


There  are  two  correlation  states.  One  where  no  extraneous  (E)  lines  are  present  (C)  in  the  cell; 
the  other  where  E  lines  do  exist  (uncorrelated)  (C)  in  the  cell  (in  one  or  both  segments).  Evaluation 
of  (2)  for  each  cell  depends  on  the  likelihood  functions  of  the  measurements  (S,  D)  and  the  states 
[S,  D,  C] .  Assuming  the  amplitude  distributions  of  target  and  extraneous  lines  are  the  same,  then 
(Appendix  El)  the  expected  likelihood  ratio  of  (2)  for  one  of  the  cells  of  type  m  (defined  below)  is 
given  by 


P(S|CD)  P(DgJA)  P(DC|A) 
P(S|CD)j  p(DCIA)  +  P(DC|A) 

*- Asm  -*  ‘ — ‘“xm2-‘ 


(5) 


where  Asm  is  the  expected  likelihood  ratio  conditioned  on  the  CD  states  of  a  cell.  The  other  terms 
are  a  priori  probability  of  the  cell  DC  states  conditional  on  associacion/nonassociation  (the  existance/ 
nonexistance  of  the  same  target  source  in  both  segments). 


The  form  of  the  terms  of  (5)  depend  on  whether  the  cell  contains  two  lines  (detection  in  both 
cells),  one  line  or  no  lines.  Using  index  m  to  make  this  distinction,  let  m  =  2,  for  two  lines,  m  =  1, 
for  1  line  and  m  =  0  for  no  lines.  Then  (5)  is  written  as  (for  the  ith  cell) 

Am  =  Asm  Xml  +  xm2  m  =  0,1>'2  <6) 

Then  the  total  likelihood  ratio  (1)  is  given  by 

n2  nl  n0 

AS  =  TT  (As2  X21  +  x22>  TT  (Asl  X11  +  x12>  TT  *0  (7) 

where, 

n£  =  number  of  cells  with  two  lines  (m  *  2) 

nj  *  number  of  cells  with  one  line  (m  *  1) 

ng  »  number  of  cells  whh  no  lines  (m  *  0) 
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D3  LIKELIHOOD  RATIO  (Agm)  DEVELOPMENT 


Asm  “  evaluated  by  integrating  over  X  and  applying  Bayes'  Rule  (Appendix  E2)  to  obtain 
P(SQ  Dq|X)  P(SjDj|X)  P(X) 


lsm 


_  ^  P(PIC) 

[/  P(S0D0IX0)  P(X0)]  [£ P(SjDj|Xj)  P(Xj)  ]  ‘  P(D|C) 


(8) 


m 


m 


If  SQ  is  not  detected  then  P(S0DQ]X)  =  P(DQ|X);  similarly  for  Sj  (Appendix  E2).  Note  for 
cells  with  no  lines  (DQ,  Dj),  (8)  reduces  to  [Asm  =  1] . 

If  only  the  detect/no  detect  event  is  “ecorded  (amplitude  not  recorded)  as  in  a  clipped  processor, 
then  'Appendix  E2)  [Asj  =  Ag2  =  1] .  Consider  each  in  turn. 


As2  —  Double  Detection 

For  detection  state  Djj  where  detection  occurs  in  both  the  o  and  j  segment  (m  =  2)  then, 

/N0NjP(X)  ,  P(D0Dj) 

As2  f No  P<X/Nj  PW  *  P(Dq)  P(Dj) 

i -  - •  * — g2 — * 

where  NQ  and  Nj  are  normal  distributions  of  the  measured  amplitude  conditioned  on  amplitude  state  X. 
These  are  likelihood  functions;  a  function  of  X.  DQ  and  Dj  are  the  detection  events  of  segments  oj. 
P(Dq,  Dj)  for  instance  is  the  joint  probability  of  the  measured  detection  events  in  cells  o,  j  given  that  a 
detectable  target  line  exists. 


As  developed  in  Appendix  E5 


As2 


/  S0  -  Sj  \  eX  -save  -T+.5aF) 
\  \J  2  op  )  X 


(9) 


where  SQ  and  Sj  are  the  measured  amplitudes  in  segments  o  and  j,  Save  is  the  average  amplitude,  T  is 

q 

the  threshold,  op  is  the  variance  of  the  log  normal  fluctuation  distribution  and  X is  a  factor  in  the 
prior  amplitude  distribution  (E6).  It  can  be  shown  that  PF  "  Xo  where  PF  is  probability  of  fade. 


Agl  —  Single  Detections 


For  cells  with  a  detection  in  one  segment  (e.g„  o)  but  not  the  other  (j)  then; 

_  jN0P(DjlX)  P(X)  P(D0Dj)_  (1( 

Agl  f  N0P(XjP(EiiXtPfXj  P(Dq)  tH&f) 

• - Fx  - ,  u.  Gj  - 1 

where  D  is  the  measured  event  that  no  detection  occurs  in  the  jth  segment.  The  second  intergral  in  the 
denominator  is,  by  inspection,  equal  to  P  (Dj)  and  hence  cancels  as  indicated.  Following  the  same 
assumptions  and  techniques  as  developed  in  Appendix  E6. 


Note  that  the  same  expression  applies  when  only  the  jth  segment  is  detected;  simply  replace 
SQ  with  Sj. 

A0  ~  No  Detection 

By  inspection  (Appendix  E2) 


This  completes  the  development  of  the  likelihood  (Agm)  terms  of  (7).  Now  consider  the  a  priori 
terms  (Aml,  Xm2)  for  m  =  1,2. 

D4  -  PRIAR  DISTRIBUTION  TERMS  OF  D2(6)  DEVELOPED 

These  terms  (Aml ,  Xm2 )  express  the  a  priori  probability  ratios  for  the  various  detection/ 
correlation  states  (DC)  conditioned  on  the  presence  (A)  or  non  presence  (A)  of  a  common  target  in  both 
segments. 

As  developed  in  Appendix  E7,  these  ratios  are  given  in  terms  of  the  target/extraneous  line  (“noise”) 
parameters  previously  defined  (np  ng,  Pp).  The  results  are  given  in  Table  D-l. 


TABLE  D-l  -  X  TERMS 


—  Double  Detection 

X21  =  (1  -  Pp)  ^  nE  2 
x22  =  nE  (nE  +  nT^  !  nE  2 

X  j  —  Single  Detection 
PF 

XH  =  IjT  nT  1  nE 
x12  =  nE  I  nE 
where, 

PF 

nE  =  nE  +  ^  2^  nT 

No  Detection 

no  nT  (1  _  PF) 

H  x0  =  e  (Refer  to  E8) 

where, 

no  =  nC  -  <nl  +  n2> 

Thus  the  complete  spectrum  likelihood  ratio  is  given  by  (7)  where  the  terms  A 
(9)  and  (11)  and  the  X  terms  given  in  Table  D-l. 
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APPENDIX  E  -  BACKUP  ANALYSIS  FOR  APPENDIX  D( A s DEVELOPMENT) 


El -EXPANSION/DEVELOPMENT  OF  PfSDI  Al/PfSDl  A) 

For  a  given  cell,  Ag  is  given  by  (cell  subscript  (i)  is  implied). 

P(SD|A) 

AS  =  P(SD|A)  (1) 

where  S,  D  are  measured  amplitude  and  detection  events  respectively,  and  A,  A  are  the  association, 
non  association  conditions  respectively  for  the  complete  spectrum.  Recall  that  D  can  take  on  one  of 
four  “values”;  detected  in  both  segments  Dj  jj  detected  in  only  one  (Dqj,  Djq)  or  no  detections  (DOQ). 

S  D  is  a  function  of  whether  or  not  the  cell  contains  an  extraneous  line.  If  it  does  then  the  segments 
are  uncorrelated  (C)  otherwise  they  are  correlated  (C). 

Including  the  correlation  state  variable  (C)  with  the  two  possible  values  (C,  C),  and  summing  over 
the  two  values  yields 

As  =  I  P  (SDC  |  A)  /  Z  P  (SDC  |  A)  (2) 

c  c 

Applying  Bayes’  Rule  yields 

As  =  2  P(SICDA)  P(DC|A)  /  2  P(S|CDA)  P(DC|A) 
c  c 

_  P(S|CDA)  P(DCJA)  +  P(S|CDA)  P(DC|A) 

~  P(SICDA)  P(DCJA)  +  P(S|CDA)  P(DCIA)  (3) 

Note  that  the  cells  cannot  be  correlated  if  unassociated  hence  P(DC|  A)  =  0;  hence  the  first 
denominator  term  vanishs. 

Assuming  that  target  and  extraneous  lines  have  the  same  distributions,  than  the  amplitude 
densities  (SI  -)  are  a  function  only  of  the  correlation  condition;  the  association  condition  (A, A)  being 
irrelevent  so  that  P(SICDA)  -  P(S|CD),  thus 


E-l 


A 


S 


P(S|CD) 

P(S|CD) 


P(DC|A)  P(DC|A) 
P(DCIA)  P(DC|  A) 


sm 


E2  -  EVALUATION  OF  A 


8m 


‘sm 


A  P(S|CD) 
"  P(S|CD) 


P(SD|C)  /  P(Dlg) 
P(SD|C)  /  P(D|C) 


(4) 


(1) 

(2) 


The  amplitude  term  of  the  numerator  and  denominator  are  determined  by  introducing  the  state 
variable  X  to  the  joint  density  function  and  then  integrating  it  out  to  obtain  the  desired  marginal 
densities.  Thus, 


P(SD|C)  =  /P(SD|XC)  P(X|C)  dx 
X 


P(SD|C)  =  /P(SD|XC)P(X|C)dx 
X 


Substitution  into  (2)  then  yields 

/P(SD|XC)  P(X|C)dx 

A  ,  X  _  P(D1C) 

8m  /P(SD|XC)  P(X|C)dx  •  P(D|C) 
X 


(3) 


(4) 


Recall  S  =  (SQSj).  For  the  correlated  condition  (C)  a  single  state  (X)  applies.  For  C,  different  X 
apply  (XQXj)  so  that  (dropping  the  correlation  notations) 


"\m 


/ P(S0D0IX)  P(SjDjlX)  P(X)  dx 


^p(SoDo|Xo)  P<Xo>  dxo  £P<SjDj'Xj)  p(Xj)  dxj 


X. 


P(DIC) 

P(D|C) 


(5) 


Consider  the  detection  and  no  detection  condition  that  may  be  measured  for  the  o  or  j  segment. 
Applying  Bayes’  Rule  for  segment  o  (similarly  for  j)  yields  the  alternate  expressions 


p(S0D0IX) 


P(D0IS0X)  P(SQ|X) 
P(D0IX)P(S0ID0X) 
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(6) 

(7) 


If  detection  occurs  (D  =  IX)  then  applying  (6) 


P(D0IS0X)  =  P(D0IS0)  =  1 


thus, 


(8) 


P^DJX)  =  P(SGIX) 


If  no  detection  occurs  (DQ  =  DQ)  then  applying  (7) 


P(S0ID0X)  =  P(S0ID0)  =  1 


thus. 


(9) 


P(SqD0IX)  =  P(Dn|X) 


In  conclusion,  the  result  (Asm)  is  given  by  (5)  noteing  per  (8),  (9)  that 


P(S0D0IX)  =  P(SQ|X)  for  the  detection  condition 

=  P(DQ|X)  for  the  no  detection  condition 
The  same  applies  for  segment  j. 

Clipped 


If  only  the  detection  event  is  recorded  (no  amplitudes)  then 


P(S0DQ|X)  =  P(D0!X) 


(10) 


or. 


(11) 


P(SjDjlX)  =  P(DjlX) 


Then  the  numerator  integral  of  (5)  equals  P(DQDj|C)  or  simply  P(D|C)  and  the  denominator  integral 
equals  P(DQ)  P(Dj)  or  simply  P(D|C).  These  terms  then  cancel  with  corresponding  terms  so  that 


Asm  ”  * 


(12) 
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E3  DETECTION  TERM  DEVELOPMENT  WITH  Q  APPROXIMATIONS 


Several  detection  terms  will  be  developed  as  given  by 


1.  P(DQ)^  =  Probability  of  detection  given  the  existance  of  a  detectable  line  i;i  a  single 

segment  cell. 

2.  P(D0Dj)  =  P(D0Dj|C) 

=  Joint  probability  of  detection  given  the  existance  of  a  correlated  line 
(same  target)  in  a  common  cell  of  both  segments. 

3.  P(D0Dj)  =  P(DQDj|C) 

=  Joint  probability  of  detection  in  one  cell  and  no  detection  (D)  in  the  other 
given  the  same  conditions  as  in  2. 


P(D)  is  obtained  by  introducing  the  state  variable  (X)  and  integrating  it  out  so  that 

P(D0>  =/p(D0X)dx=/p(D0IX)P(X)dx 
X  X 

Similarly, 

P(D0Dj)  =  /p(DqIX)  P(DjiX)  P(X)  dx 
P(D0Dj)  =  }p(D0IX)P(Dj|X)P(X)dx 


(1) 

(2) 

(3) 


P(Dq|X)  and  P(DQ|X)  Development 

P(DQ|X)  is  the  probability  that  the  amplitude  (SQ)  exceeds  a  threshold  (T),  conditioned  on  X. 
The  conditional  probability  density  function  of  amplitude  is  P( SIX)  so  that 


P(D  |X)  =  J  P(S|X)  ds 
S  =  T 

Assuming  a  Gaussian  fading  model  with  variance  as  previously  discussed,  then 


(D  P(Dj)  -  P(DC)  and  P(DQDj)  -  P(D0D.) 
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(4) 


(5) 


Then,  substitution  into  (4)  and  defining  P(D|X)  =  Q  or  Q(X)  yields 


A  ,  /s- x 

Q  =  P(Dq|X)  =  J  N 1 -  ds 

°  t  \  a  / 


S  =  T 

Note  that  this  has  the  form  of  an  error  function  (ERF).  For  X>T,  a  good  approximation  is 

r  8  f x  ~  t  '  i 

Q  =P(D0IX)=  |_1  -  .5  e  \  a  J  J 

■«(“) 

Since  P(DQIX)  =  P(Dj|X),  then  P(DQ|X)  •  P(Dj|X)  =  Q2.  An  approximation  to  Q2  is  obtained  by 
inspection  of  plots  of  Q  and  Q2  (Refer  to  Figure  El).  Q2  is  essentially  the  same  as  Q  with  an  X 


(6) 


(7) 


shift  of  about  0.5o. 


Thus, 


Q 


/X-T-.5o\ 
a  Q\  «  ) 


Thus,  from  (7) 

Q2 


P(DD|X) 


/X-T  -  .5 a  \ 

- J 


(8) 


(9) 


Note  that  if  a  line  is  present  (detectable)  it  is  either  detected  or  not  detected  so  that 

P(D|X)  =  1-P(D|X) 

=  1-Q 

Substitution  of  (7),  (9)  and  (10)  into  (1)  -  (3)  yields 


(10) 


P(D  )  =  P(Dj)  =  /q  P(X)dx 

*  V 


P(D0Dj)  =  P(D0Dj) 


JQ2  P(X)  dx 
X 

/q(1  -  Q)  P(X)  dx 
X 


(ID 


where  Q  and  are  given  by  (7)  and  (9). 
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E4  -  EVALUATION  OF  NA 

The  expression  to  be  developed  is  given  by 

na  =  /ryMx 

x 

where  NQ,  Nj  are  normal  distributions  of  SQ,  Sj  with  true  value  X.  Thus,  for  equal  variances 


Refer  to  Figure  E2. 


Completing  the  square  yields 


Since  the  first  exponent  is  independent  of  X,  it  can  be  removed  from  the  integral  so  that 


Since  the  integral  is  a  normal  distribution,  integration  yields  unity  so  that 


For  the  more  general  case  where  the  c’s  are  unequal,  define  a  =S0,  b  =  Sj,  then 


f  /a  -  x  \  ^ 

/b  _  X 

-V4  (- - ) 

Je  \  °a  ) 

•  6  \  ^  J  dx 

X  y/~2lf  <7a 

\f~2n  ab 
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(2) 


(3) 


(4) 


completing  the  squares  yields 


N 


A 


where. 


aab 


a 


/s 

x 


thus. 


1 


The  integral  is  also  unity,  thus  (6)  reduces  to 


E5  -  AS2  DEVELOPMENT 
aS2  =  F2/G2 


¥ 2  Development 


F 


2 


/N  N  P(X)dx 
X 


^N0P(X)dx  P(X)  dx 


Consider  the  numerator  of  (1)  and  Figure  E2. 


E-10 


where 


X 


(14) 


Substituting  NA  from  E4(4)  yields 


„S2 

V  2»r  V  2  o  X 


E6  -AS1  DEVELOPMENT 


-  Fi/G, 


Fj  Development 


/n  (1  -  Q-)  P(X)  dx 

V  J 


rN0P(S0)dx 


As  shown  in  the  development  of  F2,  the  denominator  is  P(SQ).  Assuming  that  P(X)  is  relatively  flat, 
then  PrfX)/P(SQ)  =  1  so  that,  (refer  to  Figure  E3(a)) 

F1  “  f"0  U  -  Qi>  dx 

X 


=  i-/n0Qj 


Note  that  Qj  is  the  integral  of  a  normal  distribution  (Nj)  integrated  from  the  threshold  (T)  to  ' 
(probability  of  detection  given  X).  Then,  substituting  this  integral  and  reversing  the  order  of 
integration  yields 


Ff  =  1 


dx  ]  dS: 


As  shown  in  E4  the  inner  integral  [  J  is  a  normal  distribution  of  the  amplitude  difference 


(Sj  -  SQ)  so  that 


so"Sj  \ 

“TtT  dsi 
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i 


The  integral  is  a  Q  function  so  that 


Fi  =i-Q-^ 


Then,  using  the  Q  approximation  of  E3  yields 


-  0.5  e  °'8  (727) 


Gj  Development 

„  P(DD) 

1  P(D) 

Using  the  detection  terms  of  E3  yields 


f  Q(l-Q)  P(X)  dx 

Q  _  _ 

1  ~J" QP(X)  dx 


J~(l  -  .5  e'-®  *)  P(X)  dx 


where 


=  (X-T) 


Assume  that  P(X)  is  essentially  linear  in  X  and  note  that  Q  (1  -  Q)  is  symmetrical  about  X-T 
(refer  to  Figure  E3(b)).  Then  P(X)  can  be  replaced  by  P(T)  and  integration  yields  the  numerator  term 


Gj  (num)  *  2  P(T) 


-.8  x/o 


■1.6  x/o 


.76 

-  —  a  P(T) 
a  0  P(T) 
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Consider  the  denominator  term  of  (8).  As  developed  for  G2  using  the  broad  P(X),  Q  step 

XX 

function  approximation  and  P(X)  =  Xe’A  ,  then. 


G1  (den)  =  /  P(X)  dx 

=  e‘XT  =  P(T)/X 

Finally,  since  Gj  *  Gj  (num)/Gj  (den),  then 


(ID 


Gj  —  Xo 


(12) 


E7~  A  PRIORI  PROBABILITIES  -  X 

The  problem  here  is  to  determine  the  a  priori  probability  ratios  (Xmj,  Xm2),  m  =  0,1,2.  These  X 
are  a  function  of  the  probabilities  P(DQDj  C|  A)  and  P(DQDj  C|A). 

Consider  first  the  association  condition  (A)  where  a  target  and  hence  correlated  lines  can  exist  in 
the  cells.  The  approach  is  to  develop  a  probability  tree  enumerating  all  target  and  extraneous  line 
conditions  and  the  probabilities  of  each. 


Define  T,  T  as  the  existance  (detection)  of  a  detectable,  not  detectable  target  line  respectively. 
Similarly,  for  extraneous  lines  (E  E ). 


If  EQ  is  dominant  over  TQ  (extraneous  amplitude  exceeds  target  amplitude),  then  TQ  is  treated  as 
a  non  existing  line  (TQ).  Probability  of  this  occurrence  is  0.5  assuming  T  and  E  have  the  same 
amplitude  distributions. 


The  probability  tree  enumerating  all  state/substates  and  their  transition  probabilities  is  shown 
in  Figure  E4. 


Definitions 

n,p  >  expected  number  of  detectable  target  lines 
nc  =  number  of  cells 


E-15 


rig  =  expected  number  of  extraneous  lines 
Pp  =  probability  of  fade 

At  the  initial  node,  as  indicated,  n^,/nc  is  the  probability  of  a  target,  (correlated  line)  being 
detected  in  either  segment.  Next,  (1  -  Pp)  is  the  probability  of  no  fade  given  a  line  and  hence  detected 
in  both  segments  (TT).  Given  this  condition,  the  probability  of  being  in  each  of  the  four  E  states  is 
indicated.  Note  the  probability  of  no  extraneous  lines  (EQEj)  is  ^1  - _ !i  j  .  For  the  other  three  E 

states  (an  extraneous  detection  in  one  or  both),  an  additional  branch  is  shown.  If  the  detectable 
extraneous  line  exceeds  threshold  (E  >  T),  it  is  assumed  that  it  will  mask  the  correlated  target  line 
thereby  resulting  in  the  detections  being  uncorrelated  (DDC).  As  indicated,  the  probability  of 
detection  given  a  detectable  E  line  is  .5.  The  other  branches  and  probabilities  were  similarly  developed. 

The  final  states  are  then  grouped  as  indicated.  The  first  group,  for  example,  applies  to  a  correlated 
detection  in  both  cells  (DDC). 

The  probability  of  occurrence  of  each  of  these  final  (detection/correlation)  states  is  obtained  by 
adding  the  probability  of  each  path  indicated,  where  the  path  probability  is  the  product  of  all  the 
branch  probabilities  back  to  the  start. 

The  resulting  probabilities  are  given  in  Table  E-l.  Recall  that  all  of  the  above  applies  to  the 
association  condition,  i.e.,  the  existance  of  a  correlated  target  in  the  two  segments. 

Nonassociation  Condition 

The  comparable  probabilities  under  the  nonassociation  condition  (A)  are  developed  by  a  straight* 
forward  modification  of  the  previous  results  for  the  A  condition. 

In  each  segment  the  expected  number  of  target  lines  is  given  by 
nTo  =  nT  ^  "  PF/2) 


E-l  7 


TABLE  E  l  -  DDC  A  PRIORI  PROBABILITY  EXPRESSIONS 
FOR  ASSOCIATED  CONDITION 


P(DDCJA) 


-a-pF> 

c 


"  E  \  2  nE  /  nE\  /  n  e\  ' 

+  2(.5)—  1—  +  (.5)(.5)l  — 

c  /  nc\nc/  \nt/ 


/  np  \  2 

(1-PF)  1-.5— 
nc  \  nc 


(i  -  pF> 


P(DDC|A)  = 


nE 

n„ 


I 

1  n„ 


1 - (1.25  -  .25  Pp) 

nc 


nt 


a  (nE + nT) 


P(DDCIA)  = 


nj 

_pp 

nE  I 

1-1.5 -  +  .5 

l21 

nc 

2  1 

nc  ' 

\n<) 

1 

fjrp 

n„ 


F 

2 


P(DDC|A)  = 


nT  f  PF  \lnE  (  nE  \ 

v-t)  -zV-vJ 


P(DDCIA) 

P(DDC|A) 

P(DD|A) 


=  P(DDCI  A) 
=  P(DDCIA) 
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Add  this  to  ng  to  obtain  the  new  ng  as  given  by 


A  , 

nEo  =  nE  +  nT  ~  PF/2^  =  nE 


and  set  n^  =  0  in  the  previous  association  analysis.  The  resulting  probabilities  are  given  below  in  Table  E-2. 


TABLE  E-2  -  DDC  EXPRESSIONS  FOR  N  ON  ASSOCIATED  CONDITION 


P(DDC|A)  =  0 


/ng'\2 

P(DDC|A)  =  /  —  J 


P(DDC|A)  =  0 


-  ng  /  nT 

P(DDCIA)  =  -  1-  — 


-  /  ng  \ 2 

P(DDCIA)  =  1  - - 


As  previously  defined,  the  A’s  are  defined  as  follows  and  the  above  probabilities  substituted  to  obtain 
the  results  of  Table  E-3. 


TABLE  E-3  -  X 


P(DDCIA) 

P(DDCIA) 


nT  ~  Pp)  nr. 


P(DDCIA)  nE  (nE  *  nT> 
P(DDC|A)  2  (ng')2 

PQDQA)  _  nTPF 
P(DDCIA)  =  2nv' 
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TABLE  E-3  -  X  (CONT’D) 


P(DDC|A)  =  ^E_ 

Xl2  ”  P(DDC|A)  =  nE' 

P(DD|A)  nT  t1  "  PF> 

X°  P(DD|A)  1+  nc 

no 

E8- APPROXIMATION  TO  II  \Q 

The  problem  is  to  obtain  an  approximation  for  F  =  [|  A0.  As  given  by  the  final  expression  of  E7 


\Q  =  1  +  x 

where 

nrj.  (1  -  Pp) 

x  = - 

nc 

thus 

F  =  TT°X0  =  ^on° 

=  (1  +  x)no 

where  nQ  is  the  number  of  cells  with  no  detections  as  given  by 

no  =nc-<nl-n2) 

Taking  natural  logs  (In)  yields 


(1) 


(2) 


(3) 


(4) 


InF  *  nQln  (1  +  x)  (5) 

Expanding  In  (1  +  x)  in  a  Taylor  series  noting  that  x  is  small  (of  the  order  of  0.01)  yields 

InF  =  nQ  [x  -  x2/2  +  ...]  (6) 


E-20 


o 

Neglecting  x  and  succeeding  terms  and  substituting  (2)  and  (4)  yields 


InF  =  [nc  -  (nj  +  n2)l 
-  nT  {1  -  Pp) 


nT  -  Pw) 


Taking  the  antilog  yields 


F  =  enT  (1  “  PF) 


E9 -EVALUATION  OF  Ac, 

_ _ “tot 

The  total  frequency  likelihood  ratio  for  cells  which  match  (detection  in  both  segments)  is  given  by 


As2tot=  AS2*21  +x22 


aS2  *s  given  by  (15)  of  E5  and  k21*  ^22  ^  Table  E3  of  E7.  Substitution  yields 


AS2tot=  AS2 


(1  -  Pp)  n,fnc 


nE  (nE  +  nT^ 


[nE  +  (1  -  Pp/2)  nT) 4  I  [nE  +  (1  -  Pp/2)  nT] 2 


For  insight  assume  that  Pp  is  small  and  define  rc  and  rE  as  the  normalized  values  with  respect  to  nn 


*VnT 


ng/nj 


VS2  r  +  7T — 

(l+rE)2  1  +  rE 
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Clipper  Processor 


For  the  clipper  processor,  Ag£  =  1.  Let  rc  =  200.  Note  that  X22  is  significant  only  for  large 
values  of  tE  and  then  it’s  value  is  unity.  Thus, 


kS2 


tot  (1  +rE)2 


+  1 


(4) 


ncn,p 


+  1 


(nT  +  nE)^ 

Amplitude  Processor 

If  amplitude  is  available,  then  as  given  by  (15)  of  E5,  the  match  likelihood  ratio  is 


where. 


'S2 


AS 


6  S' 


e'1/,As2  ePp  (5  S  +  .5) 


sj  Z  Tt  y/~2  Pp 


sQ-si 

— = -  =  normalized  amplitude  difference 

v  2  op 


(5) 


So  +  Sj 


-T 


op 


=  normalized  average  amplitude  relative  to  threshold 


Note  the  identity  substitution  Pp  =  Xop  has  been  made. 
Then, 


kS2 


tot 


S2  o 

(l+rE)2 


+  1 


(6) 


Figure  E5  is  a  plot  of  A52  vs  Pp  f°r  various  AS,  6  S.  Note  that  for  practical  fade  probabilities 
say  0.1  and  0.3,  A  is  fairly  sensitive  to  Pp. 

Figure  E6  shows  sensitivity  to  AS  and  various  5 S  with  Pp  set  at  0.2. 


Figure  E5  —  Ag£  vs  Pp 


Figure  E6  —  Ag2  vs  AS 
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APPENDIX  F  -  LIKELIHOOD  RATIO  FOR  BINARY  ASSOCIATION  - 


HYPOTHESIS  TEST  WITH  LINE  DETECTIONS  AND  MATCHES  AS  MEASUREMENTS 


—  An  Initial,  Simplified,  Frequency  Likelihood  Approach  — 

INTRODUCTION 

In  this  appendix  there  is  presented  a  combinatorially  oriented  derivation  of  the  likelihood  ratio 
for  testing  the  association  between  track  segments  based  solely  on  measurements  conveying  the  detection 
or  nondetection  of  matching  tonals. 

The  treatment  differs  most  importantly  from  that  found  in  Appendices  D  and  E  in  that, 

i.  The  a  priori  information  is  given  a  different  parametric  description, 

ii.  extraneous  lines  are  here  not  permitted  (false-detection  probability  from  noise 
alone  is  assumed  zero),  and 

iii.  the  amplitudes  of  the  detections  do  not  here  constitute  part  of  the  measurements. 

As  a  consequence  of  iii,  the  development  here  should  be  compared  to  the  “clipped  case”  development 
of  Appendices  D  and  E.  In  making  that  comparison,  it  is  seen  that  in  both  cases  the  log  likelihood  ratio 
is  a  linear  function  of  the  same  two  measurements;  the  number  of  matched  lines  between  segments,  and 
the  total  number  of  lines,  matched  and  unmatched,  in  the  two  segments. 

MODEL 

A  source  is  characterized  by  ng  tonals.  The  receiver  has  nc  FFT  bins;  at  times  ty  and  (later)  tz,  ny 
and  (later)  nz  tones  are  detected,  and  n2  of  the  ny  and  nz  “match”  in  frequency.  Hypotheses  are: 

Hj  —  Same  source  at  ty  and  tz 
HQ  —  Different  sources  at  ty  and  tz. 

Note  in  the  previous  analysis,  tQ  and  tj  were  used  rather  than  tz  and  ty,  respectively. 
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Likelihood  ratio  is 


PH1  ("2»ny»nz> 
A(n2’ny’nz)  =  PHo  (n2,ny,nz) 


(1) 


ASSUMPTIONS 

1.  ng  has  Poisson  distribution  with  mean  m; 

2.  Source  tonals  are  distributed  at  random  across  FFT  band; 

3.  p  =  probability  of  detection,  given  that  a  source  tone  is  present,  is  same  for  every  tone; 

4.  Doppler  shifts  ignored,  so  that  “match”  means  detection  in  same  FFT  bin  at  both  ty  and  tz; 

5.  Negligible  false-detection  probability  from  noise  alone;  and 

6.  nc  »  ny,  nz,  or  n2. 

ANALYSIS 

Under  Hj,  ny  and  nz  are  the  number  of  successes  in  two  independent  experiments  of  ng  Bernoulli 
trials  each,  with  single-trial  probability  of  success  given  by 

p  =  P"*  (detected  )  (2) 

for  each  of  the  ng  cells. 


Under  HQ,  ny  is  the  number  of  successes  in  nc  trials  with  single-trial  probability  of  success  given  by 


and  similarly  for  nz,  with 


where  ngy  and  ngz  each  has  a  Poisson  distribution  with  mean  m. 


(3) 


(4) 
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Note  that  under  Hj,  n2  is  distributed  as  the  number  of  red  balls  obtained  when  nz  balls  are  drawn 

without  replacement  from  an  urn  containing  ng  balls  of  which  n  are  red,  so  that 

y  i 


(5) 


while  under  HQ,  n2  is  distributed  as  in  (5)  but  with  nc  replaced  by  n 

S  *  < 


UNDER  H, 


f’Hj  <n2>  ny  nz>  ~  ^  ?  0  P  (n2,  ny,  n2,  ns) 

oo 

n .  =  °  P  ^n2,ny*  nz  *V  P(ny|ns)  P(nz|ns)  P(ns) 


(M  (ns-ny\ 

=  l  (**)  ny  n-nv/ns\  n  n  -  B  n 

ns  =  0  /  n8  V  \ny;Pyq  y\nJ  p  zq  »  2  mj 

\nz)  n 


(6) 


where 


q  =  1  -  p. 


Expanding  (6)  with 


H  =  — 

\  r  I  r!  Ik-. 


r!  (k-r)! 


and  using 


(7) 


2  - M _ fJB021<ny  +  n«-n2)emqa 

ns=°  ^ns  ~  ny  +  nz  “  n2^' 


yields 


nv  +  n 

pi  y  z 

q 


(8)  ’ 


<"<2» 


(ny  +  nz"n2)  m(q2-l) 


(9) 
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UNDER  H. 


^H0(n2*ny*nz)  =  n  ^  5  o  **  ^n2,ny,nz>  nsy’nsz) 


11  =q  P  (n2|ny,nz,H0)  P  (ny|ngy)  P  (n^)  P  (n  >  P  (nK) 
sy’  sz 


=,  £  1 
nsy,nsz  =  ® 


(ny\  (nc~ny\  nv  -  Dn  n 

\n2/  \nz  ~  n2/  <Pngy >  e  P  ^  (Pnsz> 

(i) 


7  -pn0,  n  -m  n  -m 

e  8Z  m  sye  m  8Ze 

i  n  1  n  ! 

z*  nsy*  sz* 


where  we  have  used  the  Poisson  approximation  to  the  Bernoulli  distribution  for  P  (nln0„)  and  P  (n_|n„_), 

y  ay  z  sz 


'W  =(:;)  ft)  "y  (>-?) 


nsy  \  nc  ~  ny 


/  >nv  pnsv 
(pnsv)  y  e  sy 


Evaluating  (11),  we  obtain 


(nc  -  ny)!  (nc  -  nz)-  Pny  4  nz  e2m(e  P  4> 

PHo(n2.ny.nz)  "  n2!(ny-n2)!  (nz-n2)!  (nc-[ny +  nz-n2l)!nc!  I'ny  ^  "nz^) 


where 


and  vr(v )  is  the  rth  absolute  moment  of  a  Poisson  variate  with  mean  n 


-r<P)  *  l 

k  -  0 
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so  that^*) 


Combining  (9,13-15)  gives 

(mq)ny  +  °z  em[2(l  -  e'P)  -  p(2  -  p)] 
A  (n2,ny,nz)  =  T •  _  .  _ 

(mq2)n2  (me  P)vn  (mep) 

y  nz 

where 

r  =  nc!  (nc-[ny  +  nz-n2])! 

(nc  -  ny)!  (nc  -  nz)! 

Using  Stirling’s  approximation  and 


(16) 


(17) 


(1  +  l/n)n  ~  e  for  n  »  i 
we  have,  for  nc  »  ny,  nz  or  n2, 

nc\n2 


■  -(5-) 


,  e  =  2.71828.. 


(18) 


and  (16)  becomes 


n„  »  n2 


**W*(Si) 


n„  ♦  n_ 
(mq)  y  z 


v  (me'p)  vn  (me‘P) 
y  z 


•  e 


m  [2  (1  -  e’P)  -  p  (2  -  p)] 


(19) 


(*> 


Burington  &  May,  “Handbook  of  Prob.  &  Stat.”  Handbook  Pub.,  Sandusky  Oh  1953- 
Pg.  77. 
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